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The impurity-induced quasiparticle interference(QPI) in the parent compounds of iron-pnictide superconduc- 
tors is investigated based on a phenomenological two-orbital four-band model and T-matrix method. We find 
the QPI is sensitive to the value of the magnetic order which may vary from one compound to another. For 
small value of the magnetic order, the pattern of oscillation in the local density of states (LDOS) induced by 
the QPI exhibits two-dimensional characteristics, consistent with the standing wave state observed in the 1111 
compound. For larger value of the magnetic order, the main feature of the spatial modulation of the LDOS is 
the existence of one-dimensional stripe structure which is in agreement with the nematic structure in the parent 
compound of the 122 system. In both cases the system shows Ci symmetry and only in the larger magnetic 
order case, there exist in-gap bound states. The corresponding QPI in g-space is also presented. The patterns 
of modulation in the LDOS at nonzero energies are attributed to the interplay between the underlying band 
structure and Fermi surfaces. 

PACS numbers: 74.70.Xa, 75.10.Lp,75.30.Fv 



I. INTRODUCTION 

One of the intriguing issues in condensed matter physics 
is the recent discovery of the new families of iron-based 
high-Tf. superconductors^ Like cuprates, superconductiv- 
ity arises from electron- or hole- doping of their antiferro- 
magnetic(AFM) parent compounds. Different from cuprates 
whose parent compounds are Mott insulators, the parent com- 
pounds of the iron-pnictides are bad metals whose resistivity 
is several orders of magnitude larger than that of normal met- 
als. When lowing the temperature, accompanied by the for- 
mation of spin-density-wave(SDW) state, 7 " 10 the parent com- 
pounds of 122 (AFe2As2) undergo a tetragonal to orthorhom- 
bic structural transition while for 1111 (RFeAsO x F y ) the tem- 
perature for the structural transition is higher than that for the 
AFM transition. Neutron-Diffraction measurement*^ show 
that the magnetic moment (0.87^) per Fe at 5K in BaFe2As2 
is substantially larger than the moment (0.36^b) per Fe in 
LaFeAsO^. 

Experiments on 1111 compound^Ii show an extra larger 
hole pocket around the F point, different from other pnictides 
due to surface polarity or surface-driven electronic structure. 
It is believed that the SDW picture established in the 122 com- 
pound could apply to the 1111 compound as well. The Fermi 
surfaces(FSs) of the iron-pnictides have disconnected sheets, 
the strong nesting between the hole FSs around the F point 
and the electron ones around the M point induces the SDW 
instability. The formation of SDW will partially gap the FSs 
and will have great effect on the QPI as well as other physical 
properties. Before unveiling the origin of superconductivity, 
understanding the ordered parent compound is an important 



task. 

Scanning tunneling microscopy (STM) on pnictides has 
provided much useful information on the electronic proper- 
ties.—^ 7 - For underdoped 122 compound Ca(Fei_ x Co)2As2, a 
remarkable nematic electronic structure^ has been observed 
by means of spectroscopic-imaging (SI) STM, where a one- 
dimensional (ID) structure aligns along the crystal a-axis with 
antiparallel spins, and the QPI imaging disperses predomi- 
nantly along the b-axis of the material. While for the par- 
ent compound of 1111 system LaOFeAs, there are two types 
of surface after cleavage. STM studies have revealed that a 
two-dimensional (2D) strong standing wave pattern^ induced 
by the QPI appears in one type of surface. And the corre- 
sponding dispersions along q x and q v are similar. Diversified 
electronic structure states have been observed in various sys- 
tems^Si^i whether they play a key role in the mechanism of 
superconductivity^^ 3 , has triggered much attention and moti- 
vated our investigation. 

Previously, based on some effective models, the QPI has 
been calculated by using different methodsi2i^ In order to 
well explain the above mentioned experiments, we solve ex- 
actly the QPI patterns as well as their Fourier component 
based on a phenomenological model 3 - , which considers the 
asymmetry of the As atoms above and below the Fe-Fe plane. 
Due to the surface effect during cleavage, the heights of the 
As atoms above and below the Fe-Fe plane may not be equal 
to each other, thus this model is suitable to study the proper- 
ties of the surface layers in the iron-pnictides and should be 
more appropriate to describe the STM experiments which are 
surface-sensitive. Our study will be within the framework of 
T-matrix method adopted by previous works The most 



interesting result is that we can obtain ID and 2D modula- 
tions of the LDOS by employing the same model. The value 
of the magnetic order has great effect on the spatial modula- 
tion of the LDOS. Small value of magnetic order leads to 2D 
pattern while larger magnetic order results in ID structure. In 
both cases, the patterns exhibit C 2 symmetry. To the best of 
our knowledge, there still lack works concerning both of the 
interference patterns. 

The paper is organized as follows. In Sec. II, we introduce 
the model and work out the formalism. In Sec. EH, we show 
the property of the normal state in the presence of impurity. 
In Sec. IV, we study the modulation of the LDOS induced 
by QPI with small value of magnetic order. In Sec. V, we 
investigate the QPI in the case of larger value of magnetic 
order. Finally, we give a summary . 



H. MODEL AND FORMALISM 

We start with the two-orbital four-band tight-binding phe- 
nomenological model- - which takes the rf A . and d vz orbitals 
of Fe ions into account and each unit cell accommodates two 
inequivalent Fe ions. The reason we adopt this model is that 
such a minimal model reproduces qualitatively the evolution 
of the FSs observed by angle resolved photoemission spec- 
troscopy (ARPES) experiments 3 - 4 - - — in both the electron- and 
hole- doped iron-pnictides. Based on this model, the ob- 
tained phase diagram 39 , dynamic spin susceptibility 4 ^, An- 
dreev bound state inside the vortex core 4 *!- as well as domain 
wall structure 42 are all consistent with the ARPES 43 , neutron 
scattering 4 ^ and STM— experiments. The tight-binding part 
of the Hamiltonian can be expressed as 



^0 — 2/v jvcr(tivj 



+ h.c.) - t 2i 



■ C iva- C " 



(1) 



where i, j are the site indices, v, v = 0, 1 are the orbital indices, 
and fo is the chemical potential. t\ represents the nearest- 
neighbor (nn) hopping between the same orbitals on Fe ions, 
t 2 and f3 denote next-nearest-neighbor (nnn) hoppings be- 
tween the same orbitals mediated by the up and down As ions, 
respectively. t\ is the nnn hopping between different orbitals. 
In this paper, we adopted the hopping parameters as in Ref.l30l 
i.e. t x = 0.5ev, t 2 = 0.4?!, f 3 = -2.0fi, and f 4 = 0.04?!. In the 
following the energy and scattering potential are measured in 
unit of 1 1 . The distance between the nnn Fe ions is a and set 
as unit, which is shown in Fig. [Ha). 

With the formation of SDW, the first Brillouin zone (BZ) 
needs to be folded into the magnetic Brillouin zone (MBZ). 
The Fourier transformation of the electron destruction opera- 
tor can be written as 



J(k+Q t )R n 



1 \ 1 r ik-Ri . 



e ikR > + r l J(k+Q t )Rn (2 \ 



where N is the number of unit cells, ^ can be sublattice 
A or B corresponding to the two inequivalent Fe ions, k 
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FIG. 1 : (color online) (a) Schematic lattice structure of the Fe layer 
in the SDW state. A and B are the two inequivalent Fe ions, (b) The 
first Brillouin zone (red line) and the magnetic Brillouin zone (inner 
square). The I and III quadrants are embraced by the blue lines while 
the II and IV quadrants are embraced by the green ones. All the 
dashed lines are not included in those quadrants. The corresponding 
Qk of moving k within I to that within III' is (-n, -n). In the area 
inside the blue (green) lines the corresponding V t = 1 (-1). 



is restricted in the MBZ and Qk is chosen to be ±(tt, jt) or 
+(n, -n), depending on which quadrant k belongs to, such 
that k + Qk is in the first BZ. We have shown an example 



in Fig. Hlb). Here c° 
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where a 



-2t 2 cos k v - 2?3 cos k x , a 2 



-2?3 cos k y 

k x +k y 
~2 



2t 2 COS k x , «3 = -2t4(cosky + cos£: v ), «4 = -2fi(cos + 

k — k 1 

cos - ^2 — — ). For M k the corresponding k changes into k + Q k . 
After diagonalizing the above Hamiltonian, we obtain Ho — 
Y*iivak ^yik^lvW^yik), with the energy band indices a,/j, v 
being or 1 . The analytical expressions for the eight energy 
bands can be written as e® v (k) = \{a\ +a 2 )+(-l) v aj+(-iyr- 

to, e* v (k) = $ v (k + Q k ) with F = ^r 2 + a% Y\ = \(a x - a 2 \ 
and they do not depend on the spin. At half filling fo = -0.622 
and the canonical transformation matrix reads 
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The matrix elements are functions of k, b Q0 = X\(T + Tj), 
/>qj = X2CI4, b° w = x 2 (T l - F), b° u = xia 4 and x h2 = [2(Fi + 
r) 2 +2a 2 ^y l l 2 are the renormalization factors. Substituting Eq. 
<j4j into Eq. (0, we have 



-Avka 
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In our coordinates, the origin is located at the A sublattice, the 
configuration of magnetic order is shown in Fig. QI a )- Note 
that magnetic order in YiiA n iiA ~ n ml> varies^ with site, in = 
me 'Qk-RA f or su blattice A and me l Qf( R B+m-ym f or su blattice B. 
The experimentally observed SDW 39 term is introduced as, 
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Vj = 1 for the first and third quadrants, otherwise V* = -1, as 
shown in Fig. 1 (b). <x being + 1 corresponds to spin up and spin 
down, respectively. A single impurity V s £ vo - c[ 0vo _c A o VC r + 

V m Tiva- °~ c A0vcr CA0va ~ * s l° catec l at me origin in sublattice A, 
the Hamiltonian of the impurity part can be written as 

H imp = | Yj (Vs + r V ^WQtfv<^]M k ytfy*W, (9) 



,'a'k' 



where V s and V m represent the nonmagnetic part and magnetic 
part of the impurity potential, respectively. The total Hamilto- 
nian is H — Hq + Hs dw + Hi, np . In the following we will solve 
the QPI state. 

Define the two-point Green's function as 

k '-> = -mr^yAk, TW^tf, 0)>, (10) 

where T4>{t) denotes the Fourier transform of 4>{t) in Matsub- 
ara frequencies, T T is the time ordering operator and i/t(t) = 
e TH tf/e~ rH . By using the equation of motion for Green's func- 
tion and ?§p- = e TH [H, i//] e - TH we obtain 
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D^WJco) = J^b^ y (nG^'Z(k'',k'-ico„), (13) 
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where G^ik, icj„) = [ia)„ - e£ v (k)]~ l is the bare Green's func- 
tion. Since the translational invariance is broken by the impu- 
rity, the Green's function depends on two momenta k and k' . 
Solving the Green's function is the basis for calculating the 
LDOS, to this end, we introduce a 4 x 4 matrix S, 
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where / is the 2x2 unit matrix and 

DQ-lV/-|0 R alV G 

alv 1 alv 



(15) 



From Eqs. ( fT4b and (TT3T >. we finally obtain the Green's func- 
tion 

G TJy<f^^ > = A^,(vcr£)G° Vv (MwJcW(W£tt< 
+ ^(V s +vV m )D v a %, v , IT ,{k',ico n )f a l tlv ' T {k\ (16) 

a"fi" 

where A = § _1 is a matrix, the upper index apt denote the 
row of the elements and the lower ones denote the column of 
the elements. At this stage lvloJ (kf, ioj n ) is still unknown. 
Combining Eqs. ( fT3l and ( fToT ), we obtain 

D v ^,{k!,iio n ) = f^' v ' a '(k')[l-2(V s +<rV nl )f^r\(m 
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The poles of the Green's function consist of the poles of 
the bare G° and poles of D v J /i ,(k',ioj n ), the latter ones sig- 
nify the appearance of new bound states due to the im- 
purity. From Eqs. ([181(1211) we can see that the index 
v'cr' of D can be omitted. We also note the pole of 
D£p,(k? , io) n ) is related to the magnitude of the magnetic 
order since it contains A. If we consider only the diago- 
nal term of A, then G^ /ir ,(k,k;w) n ) = G^ v (&, ico n )8 kk , + 
G aflv (k,icj„)T„, alr G° a , M , v (k',icj„), which is in the form of 
Dyson's equation. In real space, the LDOS of each site is 
p(n,aj) = -±Xy<rm-r(T T cy iv(r (T)cl iv(r (0))l Note 

CA hv (T)4 jKr (0)=2 ^Wclt^O)^^ ^, (22) 

a'k' 

c Biw (r)4jO)=JeUr)C(0)/y, (23) 
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we derive the LDOS on sublattice A and B in real space, 
respectively. Throughout the paper we have Xa - e mrA '%, 
Xak = e*&- r ", r A = (m.na), q = r AB = (0.5,0.5), and 

«i,«2 are integers, i.e. coordinates of sites of sublattice A. 
The LDOS is obtained via analytic continuation ico n — > co + ij] 
with rj being a tiny positive number and related to the lifetime 
of the quasiparticle. After some calculations, LDOS in real 



4 



space can be expressed as follows 
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We can see that the interband scattering only exists in the 
same v channel. Elastic scattering of quasiparticle mixes 
states having same energy but different momenta. The inter- 
ference between the incoming and outgoing waves with mo- 
menta k and k' can give rise to modulation of the LDOS at the 
wave vector q-k-k'. Such kind of interference pattern can 
be observed in SI-STM nowadays. The Fourier component of 
the LDOS (FC-LDOS) can be written as, 



Pqk 



+ lm{~ g2 (kk')MA6 k ; q + S k t q ) + + 
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where ~ gl Qc) = gtfik) + gtf (k), g 2 (kk') = g™{k)g™ \W) + 
gTiQglfik') and 
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gtffr) = 2 A^,G^, y (^^ n )^ y (^ +lv+1 (^),(28) 
= |(V s + cry m )Xi/r V "(^;v«> (29) 

Since the QPI in our model has C2 symmetry, which will be 
seen clearly in the remaining paper, the last line of Eq.(l26b 
will be zero, thus we only show the absolute value of the 
real part in the corresponding figures of FC-LDOS. The map 
is confined in the first BZ and we perform our calculation 
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FIG. 2: (color online) Panel (a) plots the Fermi surface of the tight- 
binding model in the first BZ. Panel (b) plots the LDOS of the tight- 
binding model. Panel (c) shows the LDOS of m = 0.2, 0.6 without 
the impurity. Panel (d) shows the zero temperature spectral function 
for m = 0.2 in the MBZ (denoted by grey color). There are four 
small pockets (Pvi, Pvq, Pv' q , Pv'j) aligning along the diagonal direc- 
tion and two high-intensity squares (Lv , Lvi). Panel (e) is similar to 
panel (d), but for m = 0.6, here only four pockets are left. 



(26) 



with N = 800 x 800 unit cells. We neglect the component 
of g i(k) since we want to see the QPI induced by impurity 
clearly. Different from the superconducting phase in which 
magnetic impurity and nonmagnetic one have distinct effect 
on the LDOS-S In the SDW state our calculations show that 
the effect of a pure magnetic impurity (crV m ) is very similar 
to that of a pure nonmagnetic one (V s ). While for the mixed 
scattering potential V s + o~V,„, the effect of the magnetic part 
is similar to varying the value of magnetic order, so in the 
following, we consider only the QPI induced by nonmagnetic 
impurity with different values of magnetic order. 

The impurity effect in the SDW state depends on the detail 
of electronic structure. Since we mainly focus on the low- 
energy structures of the LDOS, the FS topology should be 
important for the results. The black lines in Figs. [2 a) show 
the FSs of the tight-binding model, where the two hole pock- 
ets centered around F point (0, 0) are associated with e\ v , and 
the two electron ones around M point n{ ± 1 , ± 1 ) are associated 
with 6qv, in which v = (1) represents the inner (outer) Fermi 
surfaces of the hole or electron pockets. Without impurity, 
LDOS is uniform and site independent, Fig. Ob) shows the 
LDOS in the normal state, different from it, two peaks show 
up in the SDW state. Both the energy gap and the height of the 
peaks are increased with the increasing of the value of mag- 
netic order which can be seen clearly in Fig.[2jc). In our cal- 
culation the quasipartical damping is 7/ = 0.01. Here we cal- 
culate the spectral function A(k, to) at to — in the SDW state 
without impurity, which is imaginary part of Green's function 
multiplied by —- and is proportional to the photoemission in- 
tensity measured by ARPES experiment. As can be seen in 
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FIG. 3: (color online)For m = 0.2, image plot of LDOS in real space FIG. 4: (color online)The FC-LDOS p q (a>) for the m = 0.2 are shown 
at selected energy for SP Vs = ± 1. The x-axis and y-axis denote the in panels at selected energy for Vs = ±1 in the first BZ. The unit of 



coordinate of real space. 
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Fig. HJd), the locations of the bright pockets align along the 
diagonal direction which are denoted by Pv^Pvc-Pv^fVj, 
respectively, and have relations with Dirac coneA In addi- 
tion, although the FSs are mostly be gapped, the gap value is 
extremely small around the F point, thus there are two high- 
intensity squares denoted by Lvo and Lv\ . While for m = 0.6, 
the two high-intensity squares around the F point disappear, 
there are only the bright spots along the diagonal direction 
and the pockets are enlarged compared to the m = 0.2 case 
which can be seen in Fig. Hfe). Actually, our starting model 
has C4 symmetry when rotating around an As ion, which will 
be broken by a single impurity on Fe atom even if there is no 
SDW order. 
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III. QUASIPARTICLE INTERFERENCE FOR SMALL 
VALUE OF MAGNETIC ORDER 



In the SDW state m = 0.2, when the SP is weak, the LDOS 
on the impurity site has finite value. Due to the scattering 
of impurity, LDOS is site dependent. Spatial modulations of 
LDOS at energies -0.047 and 0.18 are shown in Fig. [3] for 
Vs = +1. there are 33 x 33 sites in sublattice A and 32 x 32 
sites in sublattice B. Those energies corresponding to the two 
SDW peaks of LDOS for Vs = 1 at the impurity site which 
we do not show here. Intensity of LDOS is enhanced at the 
impurity site for Vs = 1 . On the contrary, it is suppressed for 
Vs = - 1 . We can see that modulations exist along x-axis as 
well as along y-axis. In the case of Vs = -l,a> = 0.18, the 
QPI exhibits a 2D ripple-like modulation, the wave length is 
about 4a. 

Then we plot the image map of the FC-LDOS in the SDW 
state at selected energies in Fig. [4] We can see that in q space, 
there appear high intensity lines, the corresponding wave vec- 
tors are responsible for the QPI. The interference is nearly 



FIG. 5: (color online)Similar to Fig[3] but for SP Vs = ±3. 



equally strong along x-axis and y-axis. At energy a> = -0.047, 
for SP Vs = +1, the high intensity lines near the center are 
due to the scattering of pockets(Pv,,Pvp to the corresponding 
squares(Lv,). As seen from Fig.|4j for Vs = — 1, a) = 0.18, the 
value of the scattering q along the circle is about \q\ = Q.5n, 
consistent with the wave length 4a in real space. It also indi- 
cates that when bias energy deviate from zero, the underlying 
band structure is very important. 

As SP increases to Vs = +3, the LDOS at the impurity 
site decreases rapidly and will vanish for larger SP. Fig. [5] 
shows that the modulation of the LDOS still exhibits C2 sym- 
metry. At u> = 0.047, 2D ripple-like pattern appears. From 
the q space map of p q {<jS) in Fig.|6j we note that for strong SP 
Vs = ±3, at the Fermi energy, the intra-pocket scattering leads 
to the two high intensity small arcs near the center. While the 
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FIG. 6: (color online)Similar to FigH but for SP Vs = ±3. 



scattering between Lvo(Lvi) and Pvo, PVJPv\,PV^) lead to 
the off-diagonal high intensity spots. For Vs = -3, to = 0.0, 
away from the center, along the diagonal direction the high 
intensity arcs arise from inter-pocket scattering(from Pv, to 
Pv';) which can be seen clearly from Fig. [6] At a> = 0.047, the 
hight intensity wave vectors along the circle are responsible 
for the ripple-like modulation in real space. The interplay of 
FS with the underlying band structure has crucial effect in the 
scattering process, since the high intensity spots are obviously 
related to the band structure. For more strong SP, the feature 
of LDOS and the modulation of LDOS is similar, we also note 
that for large value of SP, the difference between the repulsive 
and attractive potentials becomes less obvious. 



IV. QUASIPARTICLE INTERFERENCE FOR LARGER 
VALUE OF MAGNETIC ORDER 



Previous discussions show that the value of magnetic order 
has high influence on the spectral function, thus we expect it 
will affect the QPI as well. For larger magnetic order m = 0.6 
and weak SP, at the impurity site, the asymmetry of the LDOS 
is remarkable. For positive SP Vs = 0.5, the negative energy 
peak of the LDOS is much higher than the positive one. While 
for Vs = 1 the resonance peak is enhanced and pushed to 
u> = -0.063, near the Fermi energy. Fig. |7]shows it clearly. On 
the contrary, for negative SP, the intensity of the overall LDOS 
is relatively small and the positive energy peak is higher. The 
most striking feature of the larger-m system is the existence of 
ID modulation of the LDOS. Compared to the positive Vs = 
1 case, the ID structure is more remarkable for Vs = -1, 
as can be seen in Fig. [8] For Vs = -1, at selected energies 
-0.063, the ID stripe pattern is pronounced. The existence 
of ID stripe is consistent with the nematic electronic structure 
observed in the parent compound of the 122 systems 18 which 
have large magnetic moment. 

For m = 0.6, only four pockets are left in the plot of spectral 
function, which play a important role in the formation of stripe 



= 2 

IB 

^3 



(0,0) 


--- bulk 




— Vs=0.5 




VS--0.5 











-1 





CO 



(0,0) 


-- bulk 


— Vs=1 


1 


Vs=-1 




i 



FIG. 7: (color online)Panel(a) shows LDOS at the impurity site (0, 0) 
for different scattering potential Vs = ±0.5 with m = 0.6. Black 
dashed line represents LDOS of the bulk, the red solid(blue dot) line 
corresponds to positive(negative) SP. Panel(b) is similar to (a) with 
Vs = +1. 
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FIG. 8: (color online)Similar to FigfJ] but with m = 0.6. 



patterns. Responding to the modulation of the LDOS in real 
space, dispersive excitation in q-space should appear along 
its perpendicular direction. This is illustrated in Fig. |9j we 
can see that the dominant high-intensity spots are distributed 
along a diagonal direction with differently detailed patterns at 
different energies. For Vs = -1, at a> = -0.063, the scattering 
vector q (+7T, ±7r). They are corresponding to the distance 
of ~ 2a strip pattern in real space. While the intra-pocket 
scattering leads to the high intensity spots around the center. 
Except for the high intensity spots around the center, the high- 
intensity spots for Vs = -1 correspond to the dark ones for 
Vs = 1, it means that for repulsive scattering inter-pockets 
scattering is weak. 

For strong SP Vs = +3, the LDOS at some sites in the vicin- 
ity of the impurity is strongly affected. Pronounced in-gap res- 
onance peaks appear as shown in Fig. [TO] For Vs = 3, one in- 
gap peak is located at the negative energy a> = -0.031. While 
there exist two very close resonance peaks for Vs = -3 at the 
positive energies co = 0.047 and 0.094, respectively. Those 
in-gap peaks reflect the formation of bound states induced by 
QPI. We show the image map of the LDOS in real space at se- 
lected energies u = 0.047 and -0.031 for Vs = +3 in Fig.fTTI 
We can see that, the ID stripe modulations of the LDOS are 
remarkable in all cases. The sites with in-gap resonance peaks 
are located along the lines y — -x ± 2. Fig. [L2l shows p q (a>) 
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FIG. 9: (color online)Similar to Fig|4] but with m = 0.6. 
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FIG. 10: (color online)Panel(a) shows LDOS at the different sites for 
Vs = 3 with m = 0.6. Black dashed line represents LDOS in bulk, 
the red solid line represents the site (2, 0), and the violet dash dot line 
represents impurity site. Panel(b) is similar to (a) with Vs = -3. 



for m = 0.6, Vs = +3. As can be seen, at all selected energies 
the QPI wave vectors are along the diagonal direction, though 
they form different patterns. For Vs = 3,co = -0.031, the 
width of the stripe-like pattern is apparently extended since 
the underlying band structure plays a important role at ener- 
gies away from zero. For Vs = -3, a> = -0.031, the dominate 
scattering q is still about (+0.77T, +0.77r) thus the distance be- 
tween two stripes in real space is about 2a. 

For unitary case, LDOS is identity for positive and negative 
SP. Our calculations show that at low energy u> = 0.016, the 
spatial modulation of the LDOS still has stripe pattern. How- 
ever, at higher energy a> = 0.187, beside the ID high-intensity 
stripes there are two circles around the impurity. The disper- 
sion pattern evolves with energy, at high energy a> = 0.187 
the pattern changes dramatically and has 2D characteristics. 
Therefore, the dispersion relation is complex in the multi-band 
system, and can not be fitted by a simple function as in the 
cupratesii 



V. SUMMARY 

We have investigated by the T-matrix method the modula- 
tion of LDOS and FC-LDOS in the SDW state of the iron- 
pnictides induced by QPI for different impurity strength and 
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FIG. 11: (color online)Similar to Fig[U but for SP Vs = ±3. 
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FIG. 12: (color online)Similar to Fig|9] but for SP Vs = ±3. 



well explained the two experimental works [18; s 19]. QPI is 
sensitive to the value of the magnetic order which may vary 
from one compound to another. 

For small magnetic order, beside the high intensity small 
pockets aligning along the diagonal direction, the zero-energy 
spectral function exhibits high-intensity squares around the T 
point, therefore it is easy to form 2D QPI patterns. Our cal- 
culations show that the 2D patterns of LDOS exist in real- 
and q-space for various SPs. The exact pattern varies with the 
energy and in some cases QPI induces ripple-like Friedel os- 
cillations. It is consistent with what have been observed in the 
1111 compound^ 

For larger magnetic order, the main feature of the spatial 
modulation of the LDOS is the ID structure at the energies 
lower than SDW gap. The QPI pattern in momentum space 
also supports the formation of unidirectional nanostructuresi^. 
Negative SP favors the inter-pocket scattering more than the 
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repulsive one. The LDOS on some sites in the vicinity of the 
impurity shows sharp in-gap resonance peaks since the cor- 
responding ungaped Fermi surfaces are enlarged and scatter- 
ing have more probability to induce excitation at low ener- 
gies. Our calculations show that remarkable ID stripe struc- 
ture aligns along the FM direction in real space. The reason 
is that for large-m system, zero-energy spectral function has 
four isolated pockets along the AFM direction. In addition, 
the topological analysis 14811 of a two-band model and a 5- 
band model showed that the stable ungaped Fermi pockets are 
along the AFM direction thus we expect the QPI obtained in 
those models should have ID stripe structure along the FM 
direction in real space, similar to our results. 

Our model has C4 symmetry around As iron, both the impu- 
rity and SDW could reduce the symmetry to Ci. The ungaped 
Fermi pockets as well as the underlying band structure have 
contribution to the QPI at bias energies away from zero. We 



obtain the ID and 2D QPI patterns observed by experiments 
based on one phenomenological model, the microscopic ori- 
gin of ID and 2D QPI patterns is the shape of spectral function 
at low energies. 
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